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We revisited the decaying dark matter (DDM) model, in which one collisionless particle decays early into two collisionless particles, 
that are potentially dark matter particles today. The effect of DDM will be manifested in the cosmic microwave background (CMB) 
and structure formation. With a systematic modification of CMB calculation tool camb, we can numerically calculated this effect, 
and compare it to observations. Further Markov Chain Monte Carlo cosmomc runnings update the constraints in that model; the 
free streaming length zips ^ O.SMpc for nonrelativistic decay, and ( ^"y" ^ ^5 5 x 10"^' for relativistic decay. 
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1. Introduction 



Decaying dark matter (DDM) in which the decay process C cPk ^ I 
happens at an early stage of the universe is natural for many hij{x) - I ^-e'^ 'j^/^j 
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and the graviton hij can be decomposed into Fourier modes h 
and 77 
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^p^models, and it is sometimes introduced as a way to adjust the 
Cl^DM relic density, because after decay the DM relic density will 
^ baively be lowered by 
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jh{k) + ihkj - ^-5ij)6ii{k)\. (3) 



One example is in |[lll, in which gravitino, overproduced by re- 
heating after inflation, will decay into the true lightest super- 
symmetric particles (LSP) axino as well as an axion. 
, However, pure gravitational constraints for that decay pro- 
cess are less understood and sometimes even simply ignored. 
Actually the model in [ 1] with their parameters should be ruled 
but 101 • Here we will present a model independent computation, 
in which the effect is calculated from the first principle, and 
• ^ can be compared directly with cosmological observation. Our 
model have both the parent particle and the daughter particles 
^ interacting very weakly, so the effect can only be manifested in 
- - gravitational effect, such as the cosmic microwave background 
(CMB) and the large scale structure formation. Our approach is 
based on a systematic modification of the CMB codes. 



2. Principle 

cmbf astjsl and cainb[4] are CMB calculation tools which 
are based on the photon line of sight integration technique, in- 
stead of solving the Boltzmann equation explicitly. They work 
in the synchronous gauge, the metric perturbation of which is 
gauged completely into the spatial 3x3 part 
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Then linearized Einstein equation gives the equations of mo- 
tion of h and 77 components 051] 
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Here overdot ' means derivative to conformal time r. is the 
peculiar velocity which is defined by (p + p)6 = ik'6Tf, and cr 
is the shear which is defined by (p + p)cr = -(kikj - ^6ij)l,'j - 
-i%kj-\5ij){T.-\6'j'l). 

With the h and 77 metric perturbation, the Boltzmann equation 
in terms of the fractional perturbation *!* is js]] 
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where q = ap is the comoving momentum and is conserved 
in expansion if the particle is collisionless, e = -\/q^ + a^nf- is 
the comoving energy, and h is the direction of the macroscopic 
flow of the fluid. f{q) is the unperturbed partition function, 
and the real partition function can be defined with fractional 

perturbation /(^, r, x', «,) - f(q)^l -H *F(^, t, x',n,)j. Usually 

f(q) is thermal distribution such as Fermi-Dirac distribution or 
Bose-Einstein distribution, but in our DDM model for daughter 
particles it is determined by the decay process. 
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After we plug into the collision term, we will find the for- 
mal photon line-of-sight integration solution to the Boltzmann 
equation. The anisotropy Aj = ~ today is given by Jstl 
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where the optical depth is k = - JJ" dT'k(T') < and the differ- 
ential optical depth is at = atigCrT, which is the common factor 
for all collision terms. Here tig is the number density of free 
electrons in coordinate space and o-j- = 6.65 x 10"^^ cm^ is the 
Thomson cross section, fi = (k ■ h), a = and Vg is the 
electron velocity. Suffix P means polarization mode and integer 
suffix labels multipoles of spherical harmonics. 

In DDM model, metric perturbations h and 77 get significant 
contribution from the daughter particles. The perturbation evo- 
lution of the daughter particle has the same series with the mas- 
sive neutrino 
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The right hand side source term of Eq. (Hat have contributions 
only from the first three perturbation modes {Sp, dp oc *Fo, 6 oc 
and cr oc ^'2). Because all particle species such as baryon, 
cold dark matter (CDM), photon, massless as well as massive 
neutrino talk to gravity, all their evolution will be modified. 

The codes also calculate the transfer functions (TF) as an in- 
termediate step in the CMB calculations, which is an indication 
of the large scale structure. The TF is by definition the nor- 
malized ratio of perturbation growth factor, from an very early 
stage to certain late stage, the normalization is taken with a very 
large scale which is out of horizon in the whole evolution 
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3. Modification 

We introduce free parameter which corresponds to the 
daughter particles' energy density today, while still keeping the 
nondecay CDM part Q.c. In this way we can treat any combi- 
nation of decaying and nondecay dark matter For convenience 
to use Eq. ([T]i, is not the whole energy density but only the 
mass contribution to energy density, namely the kinetic energy 
of the daughter particle is not included. Therefore we should 
have at least one massive daughter particle for this parameter- 
ization. Except for an extreme relativistic decay, the kinetic 



energy contribution is small and Q^/ represents the energy den- 
sity very well. We only consider the one-to-two decay process 
and introduce two mass ratios ^ and where nipi and mp2 
are separately the masses of two product particles and m,, is the 
mass of original particle. The last free parameter is the decay 
lifetime T,/. So the complete set of new parameters includes Q.d, 

Let us go through what will be modified in our DDM model, 
compared with the standard ACDM universe. First, the de- 
cay process will affect the expansion of the universe, through 
changing the equation of state. Before decay the DDM behaves 
as the CDM with a constant equation of state co - 0, while 
after decay it does not hold and the momentum and energy of 
daughter particle is subject to redshift, as what happens to mas- 
sive neutrino. Since before decay the DDM particle can be ap- 
proximated as being at rest, given the masses of the parent and 
daughter particles the initial transverse momentum px is fixed 
for a two body decay. The comoving momentum for each indi- 
vidual daughter particle is determined only by the scale factor 
a* at which the decay happens (in this paper we will always use 
a * to denote the quantity right at decay). As the scale factor 
a{T) grows the physical momentum p(t) decreases, while pre- 
serving 



q = a(T)p(T) = a*pT. 
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With this relation the energy density and pressure can be evalu- 
ated numerically for product particles in the modification, so is 
the expansion process. 

The decay is a continuous process for the set of DDM par- 
ticles, the way for our numerical study is to discretize it into 
30 ~ 60 channels, by which we can achieve 0.5% precision 
for CMB peak height. Each channel corresponds to daughter 
particles produced in a small scale factor region a* ~ a* + da* 
and has its own series. As the universe expands the chan- 
nels are gradually filled channel by channel in order. The 
perturbation evolution is described by Eq. (iTat . the only sub- 
tlety comes through the factor the unperturbed distribu- 
tion f(q) is no longer the Fermi-Dirac distribution of neutrino, 
but determined by the decay process. A number of product 
particles proportional to dQ. will be redistributed into q space 
d^q = Anq^dq - Anp^a'^da* , so up to some factor the unper- 
turbed partition function from decay is 
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and ^TTT- can be calculated 
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dXnq Td a* 
The last issue for the perturbation evolution differential equa- 
tion set is the initial condition. The initial values of all pertur- 
bation modes should naturally inherit the values before decay, 
which for DDM they are the same as CDM and all higher mul- 
tipoles vanish. So we have 

n = 5CDM, (12) 

€>\. (13) 
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Including the daughter particles' contribution to metric per- 
turbation finishes our modification. t 



4. Sample Calculation »; 

Generally speaking, because DDM induces no modification 
of the photon-electron-baryon plasma, the CMB anisotropy 
spectrum is only directly aff'ected by the Sachs-Wolfe (SW) 
effect through metric perturbation. Relativistic particles con- 
tribute more than massive particles to the peculiar velocity and_ 
the shear, because they are not suppressed by ^ or ^ factors.f 
So the SW effect will be enhanced dominantly by the relativist 
tic decay products, and the acoustic peaks will raise and mov^ 
to higher multipoles i, as what happens to a common S W tStctf 

As for the matter power spectrum, structure formation re- 
quires DM particles to condensate into clumps, in order to am- 
plify the density perturbation. Therefore the DM particles must 
be moving slowly to be gravitationally captured, just like a slow 
incident comet will be captured by the sun to form an elliptic or- 
bit but not the fast moving one with a hyperbolic orbit. If DDM 
is dominant and nondecay CDM doesn't exist or is negligible, 
after decay on small scale the free streaming effect will prevent 
structures from growing, and the TF will be much lower than 
the nondecay case. If there is still sizable CDM, the small scale 
power is not completely erased, but only get smaller 

For certain decay lifetime between nucleosynthesis and re- 
combination, there are two different approaches which may 
make the DDM model work. One is that the decay may happens 
to the dominant part of the dark matter, but the daughter particle 
has mass close to the parent particle and is very nonrelativistic, 
at least after redshift at the late recombination and structure for- 
mation epoch, so that the decay effect is minimized. The other 
is that the daughter particle is light and can be all the way rel- 
ativistic, but only a tiny part of the dark matter today comes 
from it, so that its effect is constrained. We will separately call 
them nonrelativistic scenario and relativistic scenario. A com- 
bination of the two works for sure, but their individual effects 
are primarily interesting. 

We present our sample calculations in Fig. [T] As expected, 
as we go to more DDM compoenent, larger mass hierarchy and 
longer decay lifetime, we see greater S W effect and greater free 
streaming effect. The TF is also calculated in |6]. In models 
with no CDM we see very suppressed TF on small scales, and 
the oscillation structure of which is given by baryon acoustic 
oscillation. The two variants of "Relativistic DDM" model co- 
incide nearly exactly, which implies in that limit only JIddm 
matters. 



5. Parameter Constraint 

We have further run our modified camb as calculation tool for 
the Markov Chain Monte Carlo (MCMC) program cosmomc 
10] for cosmological parameter evaluation. We are equipped 
with all current public data except for the SDSS DR8: 4 
CMB data sets which are WMAP 7 years, ACBAR, CBI and 
BOOMERANG fioll : 3 matter power spectrum data sets which 




Figure 1 : CMB anisotropy C^'^ C^^ C^^ modes and z = total matter TF for 
various DDM model. The "ACDM" model (black) has the best fit cosmologi- 
cal parameters which apply implicitly to the following models. The "Fiducial 
DDM" model (red) has = 0.226 completely replacing CDM, two mass 
ratios -jp =0.1 and -J^ = 0, and decay lifetime Tj = lO's. Then we sep- 
arately vary Clj to be half of the fiducial value while the other half is still the 
CDM (yellow), mass ratio = 0.01 (brown) and decay lifetime Tj = lO'^s 
(pink). We also show an example of very "Relativistic DDM" model (blue), 
which has = 10"^ and CDM component the same as the ACDM model, two 
mass ratios = 10"^ (so that Hddm = 1) "^p ~ decay hfetime 

Trf = 10' s. We then vary £1,, to be lO"'' (cyan) and ^ to be 10"'' (green). 



are SDSS DR7 LRG, SDSS DR4 and 2dFGRS UJD, as well 
as supernova data. We are constrained to minimal ACDM plus 
DDM model, which is a flat universe of Q^^ - 0, no hot dark 
matter - 0, standard cosmological constant w = -1, no 
tensor mode r - Q and no spectral index running «run - 0. The 
channels we used are CMB, HST, mpk, BBN, Age Tophat Prior 
and SN, which are consistent with that model and independent 
from each other In addition to the 4 decay parameters there are 
7 other parameters subject to Monte Carlo: Q/,/;^, Q.ch^, 0, Zrei, 
n,, log A, Asz- 

Although a running for more general parameter space is pos- 
sible, limited by our computing facility we are still working in 
the nonrelativistic scenario and the relativistic scenario, without 
exploring the intermediate region. First we consider the nonrel- 
ativistic decay scenario. We further constrain that all dark mat- 
ter undergoes decay, which trade Q^/ into Q^. Then we focus 
on some certain relation between two daughter particles' mass 
for further simplification. Corresponding to [ 1 1 we do the case 
that one product is massless, the other one which we focused 
on is that the two products have the same mass. The results are 
shown in Fig.|2] 

We find the allowed region can be described concisely by the 
free streaming length constraint zips ^ O.SMpc at 95% confi- 
dence level. Here the free streaming length is defined as the 
length measured today, where a DDM particle decays exactly 
at its expectation lifetime T,/, and the daughter particle travels 
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One product massless Two products identical in mass 
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Decay Time Tj Decay Time Tj 



Figure 2: Decay lifetime and mass ratio regions for nonrelativistic decay sce- 
nario, wliere tlie DM is pure DDM and one product is massless (left) or two 
products have identical masses (right). The black gridded regions which con- 
tain 95% MCMC points (red dots) correspond to ipg < 0.49Mpc (left) or 
/IpS - O.SOMpc (right), and the green gridded regions which contain 99.7% 
points correspond to ipg < 0.72Mpc (left) or zips < 0.73Mpc (right), compar- 
ing to the yellow gridded regions of previous constraint /ipg < IMpc 
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Figure 3: Decay lifetime and Hqui^ region for relativistic decay scenario. Here 
we take = 0.5. The black line which contains 95% MCMC points (red dots) 
is the constraint ( "^g™ Y)^ ^ < 5.3 X 10"', and the green line which contains 
99.7% points is (^^^yfj^ < 1-5 X 10"'', comparing to the yellow line of 
previous constraint {^^^WYf^ < 7.5 X 10"^ d. 



way for expression is OooM/i^ = 274.2^^7, where Y = '-^ 
and riooM and s are separately the DDM number and entropy 
densities. 

The allowed region can be analytically expressed as 
(^W^)^^ < 5x10-5 at 95% confidence level, which is shown 
in Fig.|3] Thanks to the high precision WMAP data, it improves 
the previous one by 3 orders of magnitude. 

6. Discussion 

We see no supporting evidence to introduce the DDM, min- 
imal ACDM is still among the best fit. We are not using the 
small scale such as the halo structure data, so we do not face the 
small scale structure problem of [8]. Here the scale of structure 
discrepancy is further constrained to be below 0.5Mpc. 

From high energy physics model building perspective, our 
results put interesting constraint on the gauge mediated super- 
symmetry breaking (GMSB) scenario in general, in which grav- 
itino is usually a very light LSP. Generically in that kind of 
model the gravitino mass is in hierarchy with all the other LSP 
candidate, and combined with large reheating temperature af- 
ter inflation it is overproduced, probably by many orders. Our 
results covers all the gravitino related decay process which is 
suppressed by the Planck scale, where gravitino is either the 
parent particle or one of the daughter particle. It is very hard to 
introduce a decay process which can solve the so called "grav- 
itino problem" and does not contradicts with observation, and 
the constraint is done in the most inevitable way, that the pro- 
cess can even escape electromagnetic and hadronic constraints. 

Going beyond the gravitino model in GMSB, this constraint 
can still judge other exotic models which satisfy collisionless 
daughter particles condition. Moreover, we would argue that 
this pure gravitational bound should apply to any DDM model 
even with interactions other than gravitation, since gravitational 
effect should always be the weakest. 



before matter-radiation equality. Normalized to ao = 1 , we have 

f""' r'"i q da 

A,s= v(T)dT= — .(14) 

In the identical products' masses case the free streaming length 
is the same for the two particles; and in the one massless prod- 
uct particle case it is defined for the massive particle, since it 
contributes to the structure formation. In fact we also checked 
the intermediate case such as one daughter particle is three 
times heavier than the other, in which the free streaming length 
cannot be defined without ambiguity, and the numerical allowed 
region is a reasonable average of the two particles' regions cor- 
responding to the same free streaming length at the same confi- 
dence level. 

Then we consider the relativistic decay scenario. Here we 
use Eq. ([TJ to trade into Qddm of parent particle if it doesn't 
decay, by doing so we are able to explore a large hierarchy re- 
gion of relativistic mass ratio, the exact value of which has little 
efi'ect as shown in the previous sample calculations. Another 
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